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Abstract. We exploit the 14 gamma-ray bursts (GRBs) with known redshifts z and the 7 GRBs for which there are constraints 
on z to determine the GRB rate ^grb( z )' usm 8 a method based on Bayesian inference. We find that, despite the qualitative 
differences between the observed GRB rate and estimates of the SFR in the universe, current data are consistent with R GRB (z) 
being proportional to the SFR. 



INTRODUCTION 

There is increasing evidence that GRBs are due to the 
collapse of massive stars (see, e.g., [1] for a discussion of 
this evidence). If GRBs are indeed related to the collapse 
of massive stars, one expects the GRB rate to be roughly 
proportional to the SFR. However, the observed redshift 
distribution of GRBs differs noticeably from that of the 
SFR: the observed GRB redshift distribution peaks at 
z « 1 and few bursts are observed beyond z ~ 1 .5, while 
the SFR peaks at z ~ 2 and 10-40% of stars are thought 
to form beyond z = 5 (see, e.g., [2,3,]). 

However, observational selection effects play an im- 
portant role in determining the observed redshift dis- 
tribution of GRBs. The important question is therefore 
whether or not the discrepancy between the observed 
GRB redshift distribution and the redshift dependence of 
the SFR is entirely due to selection effects; i.e., is the 
GRB rate roughly proportional to the SFR after taking 
observational selection effects into account? We address 
this question in this paper. 



IPN detect GRBs as a function of peak photon flux P 
by comparing the BeppoSAX and IPN peak photon flux 
distributions to that of BATSE. We fit the model jointly to 
the peak fluxes and redshifts of the 14 GRBs with known 
z, and the 7 GRBs for which there are constraints on z- 

We write the rate of GRBs that occur per unit redshift 
and luminosity as 

dN /dzdL N — p(z)f(L N ) , (1) 

where 

P(z)=R G kb(z;P,Q) x (1 +z)-' xAnr{zf{dr/dz) (2) 
is the rate of GRBs that occur at redshift z, 
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is the rate of GRBs that occur at redshift z per unit 
comoving volume (see [5]), t (z) is the elapsed time since 
the Big Bang, and 



METHOD 

We adopt a Bayesian approach. We calculate the likeli- 
hood of the data given the model, and convert it to a pos- 
terior distribution on the model parameters. We assume 
a very general model for the GRB rate, and a power- 
law model for the intrinsic GRB photon luminosity dis- 
tribution (we assume that the amplitude and the power- 
law index of the photon luminosity distribution does not 
evolve; we relax this assumption in future work). We de- 
termine the efficiency with which BeppoSAX and the 



f(L N ) = L-/ x®(L N - L mm )0(L ma x - L N ) (4) 

is the intrinsic photon luminosity distribution of GRBs. 
Thus the model has five parameters: P, Q, j3, L mia , and 

^max ■ 

We write the efficiency with which GRBs with known 
redshifts are found as s(z,P) = £ z (z) x £ ST £ p (P), where 
£ z (z) is the efficiency with which the redshifts of GRBs 
are determined from optical observations once they are 
detected by a y-ray burst instrument. We take £, (z) = 1 
for GRBs whose redshifts were determined by detection 
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FIGURE 1. Comparison of the best-fit models to the dif- 
ferential distribution of peak fluxes of BATSE, IPN and Bep- 
poSAX bursts, and the cumulative peak flux distributions for 
the three experiments. 



of an absorption-line system in the optical afterglow of 
the burst and e z (z) = ©(1 - z) + 0(z - 2.5) for GRBs 
whose redshifts were determined by measuring emission 
lines in the spectra of the host galaxy. The latter expres- 
sion accounts qualitatively for the difficulty in measur- 
ing redshifts when the H a and 0[II] emission lines from 
host galaxies do not lie in the visible spectrum. The quan- 
tity e ST is the "stereo-temporal" efficiency that accounts 
for limitations of exposure in time and solid angle and 
E P {P) is the efficiency with which BATSE, the IPN and 
BeppoSAX detect GRBs as a function of peak flux P. 
Figure 1 compare the best-fit models of e P and the cu- 
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FIGURE 2. Efficiency e(P) with which BATSE, IPN and 
BeppoSAX detect GRBs as given by the best-fit models to 
the differential distribution of peak fluxes of BATSE, IPN and 
BeppoSAX bursts. 



mulative peak flux distributions of BATSE, the IPN, and 
BeppoSAX, respectively. Figure 2 shows the best-fit e P 
for each of the three experiments. 
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FIGURE 3. Probability contours in the (P,Q)-plane for the 
GRB rate parameters (P, Q) found from fitting to the 14 GRBs 
with known z and the 7 GRBs with constraints on z. The solid 
curves correspond to the 68% and 95% probability contours. 
Also shown are the (P,Q)-values corresponding to no space 
density evolution (+), the Madau et al. SFR [4], and the Rowan- 
Robinson phenomenological model fit to IR, optical and UV 
data [5]. 

The likelihood function is then given by 

if = exp{-| dzdPn&PMz^llniz^Pi) , (5) 
where 

= p(z)xf(Anr(z) 2 (l+z) a P) x4nr(z) 2 (l+z) a 

(6) 

is the expected number of events observed within dzdP 
of (z,P). The quantity a is the burst spectral index, which 
we set equal to one in this work. By an application of 
Bayes' Theorem, we now regard if as an (unnormalized) 
probability distribution on the model parameters. 



RESULTS 

Figure 3 shows 68% and 95% probability contours for 
the GRB rate parameters P and Q. Also shown on the 
plots are the best-fit SFR models of Madau et al. [4] and 
Rowan-Robinson [5]. The SFR models lie at about a 68% 
excursion from the best-fit GRB rate model. Thus we find 
that, despite the qualitative differences that exist between 
the observed GRB rate and estimates of the SFR in the 
universe, current data are consistent with the actual GRB 
rate being approximately proportional to the SFR. 
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